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ABSTRACT 

Multi-dimensional instabilities have become an important ingredient in core-collapse supernova (CCSN) the- 
ory. Therefore, it is necessary to understand the driving mechanism of the dominant instability. Comparing 
3D CCSN simulations with turbulence theory, we find that buoyancy-driven convection dominates post-shock 
turbulence. In general, the convective fluxes and kinetic energies in the neutrino-heated region are consistent 
with expectations of buoyancy-driven convection. Specifically, the convective flux is positive where buoy- 
ancy actively drives convection, and the radial and tangential components of the kinetic energy are in rough 
equipartition (i.e. K r ~ Kg+K^). Both results are natural consequences of buoyancy-driven convection, and 
are commonly observed in simulations of convection in other contexts. Most compelling, though, is the con- 
sistency between 3D CCSN simulations and predictions of neutrino-driven convection theory. For one, global 
buoyant driving is balanced by global turbulent dissipation. Secondly, the convective luminosity and turbulent 
dissipation are linearly proportional to the driving neutrino power. Thirdly, we accurately calculate the shock 
radius only if we include turbulent ram pressure in the shock conditions. In all, these results suggest that in 
neutrino-driven explosions the multi-dimensional motions are consistent with neutrino-driven convection, and 
there is little need to invoke alternative instabilities such as the standing accretion shock instability. 
Subject headings: convection — hydrodynamics — instabilities — methods:analytical — methods: numerical 
— shock waves — supernovae: general — turbulence 



1. INTRODUCTION 

The explosive death of massive stars, in particular core- 
collapse supernovae (CCSNe), are some of the most energetic 
explosions in the Universe, and, as such, are fundamental to 
a wide range of other astrophysical phenomena. To high- 
light a few important examples, CCSNe are a major site for 
nucleosynthesis, mark the birth of neutron stars and black 
holes, and are major contributors to galactic dynamics and 
star formation. Despite their importance, understanding the 
mechanism remains an important unsolved problem. What- 
ever the mechanism, it has long been suggested that neu- 
trinos and mu lti-dimensional instabilities play major, if not 
central, roles (lEpsteinll 979t iBethe & Wilsonlll985t iBurrowsl 



1987llWilson&Mavld 



1988; Bethe 1990; Herant et al. 199 



Benz et all 1 1994 iHerant et al.l 11994 IBurrows et al 



1995 



Janka&Mulleri 11995b IBlondin et al.l l2003t iMarek & Jankal 
2007b iMurphy & Burrowsl 12008b iNordhaus et al.l 120101) . In 
this paper, we use turbulence theory to show that the domi- 
nant multi-dimensional instability is consistent with neutrino- 
driven buoyant convection. 

Multi-dimensional simulations have long suggested that 
aspherical, nonlinear instabilities play important roles in 
aiding the delayed-neutrino mechanism toward success- 
ful e xplosions. Otherwise, except for the le ast massive 
stars dKitaura et al.l l2006t IBurrows etaD 120071) . the spher- 
ical d elayed-neutrino mechanism f ails to produce explo- 
sions (ILiebendorfer et al.l 12001 bllab iRampp & Jankal 120021: 



Buras et al.l 12003b iThompson et al.l 12003b ILiebendorfer et al.l 
20051) . Even though the importance of multi-dimensionality 
is clear, which instability dominates the aspherical motions 
has been less clear. Initially, neutrino-driven convection 
was id entified as the most relevant multidimensional insta- 
bility dBurrowsl 119871: IWilson & Mavlel 11988b iBethd 11990b 

1 Princeton University, jmurphy@astro.princeton.edu, jdo- 
lence@astro.princeton.edu, burrows@astro.princeton.edu 



iBenz et al.lfl99l: IBurrows et al1fl99l Uanka & Mulie3[l99l . 
but then idealized 2D simulations discovered a new instability, 
the st anding accretion shock instability (SASI) (iBlondin et al.l 
120031) . Both must exist at some level, but it has never been 
made clear which dominates in realistic simulations. 

Investigating the importance of neutrino-driven convec- 
tion in CCSN theory has a long history. In the earliest in- 
vestigations, it was suggested that convection expands the 
shock radius, making the gain region l arger and increasing 
net neutrino heating ( Benz et all 11994 IBurrows et al.l 11995b 
Uanka & Mulleril 19951 1 199(1 7 However, none of these investi- 



gations verified that the aspherical motions are in fact due to 
neutrino-driven convection, nor did they check whether the 
turbulent ram pressure is actuall y sufficient to produce the 
expanded radius. More recently, IMurphy & Burrowsl (120081) 
considered the global conditions for explosion and found 
that turbulence reduces th e critical neutrino luminosity for 
successful explosions (se e lYamasaki & Yamadal d2006l) and 
IMurphy & Meakinl ( 20111) for theoretical discussions). Us- 
ing 3D si mulations, 



Nordhaus et al.l (120101) found similar re- 



sults, but lHanke et al I (1201 Th conclude that another multi 
dimensional instability, the SASI, might be more important 
in aiding successful explosions. 

The SASI is an instability of the standing accretion shock 
that was first discovered in idealized simulations which pur- 
posely ne^kc^edj^uttinos to suppress convective instabil- 
ities (IBlondin et al.l 120031) . These idealized 2D simula- 
tions exhibited strong up-and-down sloshing motions of the 
shock, leading to an immediate connection to the slosh- 
ing shock motions observed in more realistic 2D simula- 
tions. Consequently, most subsequent studies focused on 
the mechanism responsible for the SASI or p ostulated that 



2003b IBlondin & Mezzacappal 12006b iForiizzo et al.l 


2006; 


Marek & Janka 


2007; Schecketal. 2008; Foglizzo 


2009; 


Sato et alj|2009 


Fernandez 2010; Hanke et al. 2011). Linear 
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theory suggests that an advective-acoustic cycle is the mech- 
anism for the SASI dGuilet & Fogrizzoll2012T: IFoglizzo et ail 
2012). These analyses show that under certain conditions an 
advective-acoustic instability in addition to the buoyant in- 
stability may operate in the core collapse context. However, 
to more easily study the SASI, most analyses used idealized 
simulati ons in which buoyancy-driven instabilities were sup- 
pressed dBlondin et al.1 120031: felondin & Mezzacappal 120061: 



iSato et al.ll2009UFoglizzoll2009l) 



Studies that focus on th e SASI and include neutrino trans- 
port or heating are few (jFoglizzo et al.l 120061 IScheck et al.l 
120081 iFernandez & Thompsonl 120091). Using a toy model 
and linear theory. IFoglizzo et ak I (120061) considered the lin- 
ear growth of convective instabilities and found that advec- 
tion can sweep small-perturbation modes out of the convec- 
tively unstable region before they have time to grow to non- 
linear amplitudes. Hence, they conclude that a negative en- 
tropy gradient is not enough to drive convective instability; 
one must also consider the ratio of the advection time to the 
local buoyancy timescale (\)- For x > 3, the linear convec- 
tive insta bility succeeds, but fo r \ < 3, the SASI dominates. 
However. lFoglizzo et al.l (120061) cautioned that this analysis is 
best suited for linear growth of small perturbations, and if the 
seed perturbations are sufficiently large , convection may en- 
sue even if \ < 3. IScheck et al.l ( 120081) investigated whether 
this condition is relevant in more realistic simulations and 
found that with small initial perturbations, the SASI initially 
appeared to dominate when \ < 3- However, after ^100 ms, 
large SASI perturbations appeared to trigger convection. With 
larger, but still modest initial perturbations O(10~ 2 ivfl con- 
vection appeared to dominate at all times. Given that large 
convective perturbations in t he progenitor (B azan & Arnettl 
ll998HMeakin & Arne"ttll2007l) will provide large perturbative 
seeds, th e latter scenario is mor e likely. Based upon the linear 
analysis, IFoglizzo et alj ((2006) conclude that "advective sta- 
bilization weakens the influence of convection on the largest 
modes," but we suggest that the multi-dimensional simula- 
tions indicate otherwise. 

In many other simulations, there are hints that buoyancy- 
driven convection dominates nonlinear motions. Two- 
dimensional and three-dimensional simulations that include 
neutrinos show prominent, positively-buoyant, high-entropy 
plumes and negatively-buoyant, low-entropy plumes at late 
times (see Figure [T). Even in 2D simulations that exhibit 
large sloshing motions of the shock, outward excursions of 
the shock are accompanied by risi ng, high-entropy plumes. 
Most recently burrows et al.l (120121) have analyzed the multi- 
dimensional shock motions in 2D and 3D simulations and 
have found that the sloshing motions frequently identified 
with SASI are suppressed in 3D compared to 2D, and the 
character of the oscillations is sensitive to the driving neu- 
trino luminosity. These results are consistent with neutrino- 
driven convection as the source for the aspherical shock mo- 
tions. The correlation with neutrino luminosity is an obvi- 
ous indicator of neutrino-driven convection. The reduction 
in the large-scale sloshing modes in going from 2D to 3D is 
consiste nt with known differen ces in turbulence between 2D 
and 3D dBoffetta & Eckell2012l) . In 3D, turbulence cascades 
to smaller scales only via a constant energy cascade. Two- 
dimensional turbulence exhibits a double cascade: an enstro- 
phy cascade to smaller scales and an energy cascade to larger 
scales. This difference naturally leads to more large-scale, 

2 iv is the radial velocity 



coherent structures in 2D (see Figure [TJ. Might this be the 
source for the apparent sloshing modes in realistic 2D sim- 
ulations? Albeit circumstantial, these observations call into 
question the assumed dominance of a SASI in CCSN simula- 
tions that include neutrinos. 

Determining which instability dominates, requires a de- 
tailed analysis of the nonlinear motions and comparisons with 
theoretical predictions. Unfortunately, a nonlinear theory for 
a SASI mechanism does not yet exist, so the predictions of 
such a theory can not be falsified. On the other hand, non- 
linear theories for buoyancy-driven convection do exist, and 
can be tested against simulations that include neutrinos and a 
realistic equation of state. In this paper, we focus on the latter 
and leave the former for future work. 

To test whether nonlinear, turbulent flows of 3D CCSN 
simulations are consistent with buoyancy-driven convection, 
we compare these simulat ions with a nonlinear th eory for 
neutrino-driven convection dMurphv & M eakin 20lTJ)- in sec- 
tion |2j we describe the 3D simulations. Then in section [3] 
we highlight some of the more illuminating elements of the 
nonlinear convection theory, and in section|4] we compare the 
predictions of this theory with the properties of the 3D simula- 
tions. Finally, in section|5j we conclude that the turbulent mo- 
tions in 3D simulations are consistent with buoyancy-driven 
convection. 

2. SIMULATIONS 

The numerical results of this paper are based upon CCSN 
si mulations using CASTR O and are similar to the simulations 
of iNordhaus et al.l (120101) . CASTRO solves the hydrodynam- 
ics equations using a Godunov-type finite-volume scheme 
where the boundary fl uxes are calculated u sing an approxi- 
mate Riemann solver (lAlmgren et al.ll20lol) . Specifically, it 
evolves the conservative hydrodynamic equations: 

a, P +v-(pu) = o, (i) 

a,(pu) + V-( /9 uu) = -VP+pg, (2) 
d,(E) + X7-[u(E+P)] = pu-g + P q, (3) 



and 



where p is the mass density, u is the velocity, g is the local 
gravitational acceleration, E is pe + pu 2 /2, e is the specific 
internal energy, P is the pressure, and q is the net heating 
and cooling. For gravity, we use the Newtonian monopole 
approximation, g = -(GM/r 2 )r, and f or pressure, we u se a 
relativistic-mean-field equation of state ( IShen et al.ll998l) . As 
initial conditions for these simulations, we use the 15-M 
progenitor model of lWooslev & Weaver! (|199 5)PI 

Following the prescripti on establishe d in 

iMurphv & Burrowsl (120081) and INordhaus et all (120101 ). 
we approximate neutrino heating and cooling with local 
prescriptions, i.e. 

100km\ 2 / T v \ 2 / T 



(4) 

where L v is the luminosity of electron- or anti-electron-type 
neutrinos in units of 10 52 erg/s, T v is the temperature of the 
neutrinos (which we set to 4 MeV for all runs), T is the local 
matter temperature, and the constants are H = 1.544 x 10 20 

3 See IMurphv & Burrowsl i2008h and lHanke eTaTl GOTTI) for representa- 
tive accretion rate history curves. 



Figure 1. Entropy color maps of 2D (left) and 3D (right) CCSN simulations. Cooler colors represent lower entropies and warmer colors represent higher 
entropies. These stills represent t he flow at 250 ms after bounce for L v = 2.1 X 10 52 erg/s. The 2D simulation has a higher proportion of coherent structures, 
which turbulence theory predicts (Boffetta & Ecke 2012). Despite the differences between 2D and 3D, both show positively (high entropy) and negatively (low 
entropy ) buoyant plumes, a strong indication of neutrino-driven convection. 



and C = 1.39 9 x 10 20 . For a derivation of these constants see 
Uankal (120011) . In this paper, we consider neutrino luminosities 
ofL v = 2.1, 2.23, and2.3. 

Absorption and emission of electron- and anti-electron- 
type neutrinos is most efficient on free protons and neutrons. 
Therefore, we weight the heating and cooling terms by the 
combined mass fractions of protons and neutrons, i.e. Y p + Y„. 
Equation [4] is an approximation that is most relevant in the 
optically-thin regime. Therefore, to suppress unphysical heat- 
ing and cooling at high optical depths, we further weight 
eq. (HJi by exp(-r), where r = J npdr is an average optical 
depth of the electron- and anti-electron-type neutrinos, k is 
the neutrino opacity, and we approximate the optical depth 
with 



4 V 4MeV 



(Y n + Yp) 



10 10 g/cm 2 



dr (5) 



To simulate the range in length and time scales encountered 
in core-collapse simulations, we use CASTRO's adaptive- 
mesh-refinement (AMR) and adaptive time-stepping capabil- 
ities. We have developed an AMR strategy to simulate the 
full dynamic range of spherical collapse, while keeping the 
run-time and memory requirements as low as possible. Over- 
all, we use 6 levels of refinement, each a factor of two smaller 
than the next largest level. The largest domain of the 3D simu- 
lations is a cube with 10,000 km on a side and has a resolution 
of 32 km at the coarsest level. To adequately resolve the proto- 
neutron star (PNS) structure, the finest level has a resolution 
of ^0.5 km out to a radius of 50 km. In between, we initialize 
the refinement level (£) to maintain a roughly constant angular 
resolution of A8 <~ 0.7°, i.e. 



max mm 



log 2 64 



40 km 



,6 ,0 



(6) 



where [ J is the floor function. 

Throughout the simulation, we maintain eq. (O as the min- 
imum resolution. In addition, we set the minimum refinement 
level to 4 everywhere the entropy is greater than 5 kb /baryon 



(where kb is Boltzmann's constant). Effectively, this extends 
level 4 refinement (^2 km resolution) to include all regions 
interior to the stalled shock. As the shock expands during ex- 
plosion, level 4 refinement expands in radius requiring ever 
greater memory. To limit the storage requirements of the sim- 
ulation we impose maximum radii for each refinement level 
via 



£max = max rnin 



log 2 64 



75 km 



6k0 



(7) 



3. RELEVANT ELEMENTS OF TURBULENCE THEORY 

Using Reynolds decomposition. IMurphv & Meakinl ( 1201 11) 
developed a consistent set of equations for the evolution of 
the post-shock flow, including the effects of buoyant-driven 
convection. Here, we highlight elements of this turbulence 
theory that can be easily compared with the properties of the 
3D simulations. 

Including the effects of turbulence, the steady-state, 
spherically-averaged conservation equations for mass, mo- 
mentum, and entropy are 

V-( Po v+(pV)) = 0, (8) 
(pu ) • Vv = - VPo + pog ~ V • (pR) (9) 



and 



(p„).V. = (f ) + ^-V-<F f >. (10) 



where (•) is an average over solid angle and approximately 
one eddy turn-over time, the subscript denotes the back- 
ground flow and the prime denotes the perturbation due to 
convection, and e is the turbulent dissipation. To avoid cum- 
bersome subscripts later, we do not use for the background 
velocity. Rather, the background velocity is v and the per- 
turbed velocity is v', i.e. u = v + v'. The Reynolds-averaged 
equations are similar in form to the usual equations of hydro- 
dynamics, except these equations have three new terms that 
are associated with turbulence. The mass equation, eq. dS), 
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includes the divergence of the buoyancy flux, (p'v'}, the mo- 
mentum equation includes the divergence of Reynolds stress^, 
R = v'jv'j, and the entropy equation includes the transport of en- 
tropy by the turbulent entropy flux, F s . For low-mach-number 
flows, the buoyant flux and entropy flux can be rela ted by 
a thermodynamic derivative (iMurphy & Meakinll201 ll) . so in 
the rest of this paper, we consider only R and F s . 

The new turbulent terms requi re additional equations to 
close the system of equations. SeeyVlurphv & Meakinl d201 lb 
for the full set, but here we discuss only the equation for R, or 
more specifically, we present the specific kinetic energy (K) 
equation, where K is related to the trace of the Reynolds stress 
by K = (l/2)Tr(R). The turbulent kinetic energy equatior0is 

d{pK) /dt + v-V (pK) + (pK) V ■ v = 

+ (p'v')-g-V-(F K )-V-(F P ) (11) 
+ (P'V-v')-p e. 

On the right-hand-side of this equation, we have the work 
done by buoyancy, turbulent redistribution by the turbulent 
kinetic energy flux, the divergence of the pressure flux (Fp = 
P'v'), work done by turbulent pressure, and turbulent dissipa- 
tion. For low-mach-number flows the pressure flux and work 
done by turbulent pressure are small, so henceforth we ignore 
these terms. 

Assuming steady-state and integrating over the entire con- 
vective volume, we find a balance between global buoyant 
driving and global turbulent dissipation: 

J (p'v')-gdV = J p edV, (12) 

where we assume zero turbulent kinetic energy flux at the 
boundaries of convection. In comparing this equation with 
3D simulations, calculating the integrated buoyant driving 
term (Wb) is trivial. We merely use the simulations to cal- 
culate the integral, W/, = J (p'v')gdV. Turbulent dissipation, 
Ek = J poedV, on the other hand requires a model. We adopt 
Kolmogorov's hypothesis, in which dissipation is set by the 

largest scales^ and is of order e ~ v' 3 /£, where v' is a typ- 
ical turbulent velocity on the largest length scale, C. For- 
mally, £ is a free para meter of the model . However, simu- 
lations of stellar models ([Arnett et al.l2009l) and core collapse 
(IMurphy & Meakinll201 lh indicate that C takes on the largest 
possible value, the radial extent of the region actively driving 
convection. For this paper, we find that setting C to the size 
of the gain region satisfies global balance. In effect, C is no 
longer a free parameter, but a condition imposed by the global 
structure. If the post-shock turbulence is driven by buoyancy, 
then the 3D simulations should be consistent with eq. (fl2l) . In 
section|4] we show that the 3D simulations are indeed consis- 
tent with eq. (fT2l . 

Next, we derive analytic relationships for neutrino-driven 
convection. First, we use the entropy equation, eq. (fTOt . 
to derive a scaling for the turbulent luminosity, TqL s , where 

4 For practical purposes, we calculate R via (pwVj) /pa, which we find 

to be nearly identical to the definition given in the main text. Note, we could 
have easily defined R as pov'v'j. However, this definition obscures the behav- 
ior of the turbulent velocities with a steep density gradient. 

5 N ote that there is an erroneous extra term in eq. (8) of Murphy & Meakin 
1 201 1). We have corrected it in this reproduction of the equation. 

6 Actually, the idea that dissipation s tarts at the la rgest scales and cascades 
to sma ller scales was first proposed by Richardson ( 1922). but Kolmogorov 
119411) established the quantitative theory that we reference in this paper. 



L s = A-Kr 2 F s . In the convective region, we assume that the 
terms on the right-hand-side of the entropy equation, eq. (TTOl i. 
are of the same order. In the limit of a zero entropy gradient 
(efficient convection) or zero background velocity, the terms 
on the right-hand-side exactly balance. Even when the back- 
ground velocity or entropy gradient is nonzero, we expect the 
terms on the right-hand-side to be of the same order. Under 
this assumption, we expect convective entropy redistribution 
(V • F s ) to scale with neutrino heating (pq /T). In other words, 
T()dL s /dr ~ A-Kr 2 pq. To express this in terms of L v , we sub- 
stitute the expression for q (eq. [4]i into this expression, and 
assume that cooling is negligibly small in the heating region. 
Because the heating term is proportional to L v K,jr 2 , we arrive 
at TodL s /dr ~ L v np, where k is the opacity to neutrinos. A 
rough integration of this last scaling leads to 

T q L s ~L v t. (13) 

Therefore, if neutrino-driven convection dominates the turbu- 
lent motions, then we expect the turbulent luminosity, TqL s , 
to be linearly dependent on the driving neutrino power, i.e. 
the neutrino luminosity times the neutrino optical depth in the 
gain region, L v t. In section |H we show that 3D simulations 
are consistent with this prediction. 

For the second analytic scaling, we show that the turbu- 
lent dissipation also roughly scales linearly with the driv- 
ing neutrino power. Equation (fT2l indicates a direct corre- 
lation between the scale of the buoyant flux, (p'v'), and the 
total turbulent dissipation, = J p^edV. The buoyant flux 
can be related to the turbulent entropy flux, F s , through a 
simple thermodynamic derivative, i.e. (p'v') w r]F s , where 
r\ = (dp/ds)p/ pis evaluated at constant pressure. Hence, total 
turbulent dissipation scales with F s . Furthermore, virial ar- 
guments suggest a relationship between gravity and the ther- 
mal energy, e.g. GM/r^k/,T. Inserting these two scalings 
into eq. (fT2l i and dropping kbT) (which is dimensionless), we 
have that TqL s ~ E^. This is not a new result, but just an ap- 
proximate reiteration that buoyant driving is balanced by total 
turbulent dissipation. However, this approximate form, com- 
bined with the scaling derived for TqL s in eq. (fT3l l, leads to a 
linear scaling between the total turbulent dissipation and the 
driving neutrino power, i.e. 

Ek ~ L Ve r . (14) 

Of these two predicted scalings, the one for TqL s (eq. [T3l 
relies on fewer approximations and is more fundamental. In 
essence, the two independent predictions of buoyancy-driven 
convection theory are that the turbulent luminosity, TqL s , 
scales linearly with the driving neutrino power, L v t, and that 
global buoyant driving is balanced by turbulent dissipation. 
The scaling for (eq. H4l . is merely a recasting of these two 
statements. Because eq. (TBI is secondary, it necessarily re- 
quires more approximations. Despite these failings it is still 
informative to see how the dissipation scales with the driving 
neutrino power. So for completeness, we compare eq. (U~4t 
with 3D simulations in sectional but we are more confident in 
the derivation of eq. (fl3l , and therefore, we emphasize the lin- 
ear scaling between turbulent luminosity and neutrino power. 

An important aspect of the core collapse problem is the 
presence of the standing accretion shock, so we also consider 
how buoyancy-driven convection affects the stalled shock 
radius. Formally, the stalled shock is located where the 
upstream and the downstream profiles satisfy the Rankin- 
Hugoniot jump conditions. For zero shock velocity, the mass 
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flux, momentum flux, and energy flux conditions are 

A[pv] = 0, (15) 

A[P+pv 2 ]=0, (16) 

and 

A[e+P/p+v 2 /2]=0. (17) 

In detail, the shock position is a nontrivial solution to a bound- 
ary value problem for p, v, and P. However, with a few rea- 
sonable approximations, the shock boundary condition can be 
reduced to one expression. First, we assume steady-state and 
that the mass accretion rate (M) is constant. Secondly, we 
assume that the upstream flow is in near free-fall and essen- 
tially pressureless. Thirdly, we assume for the purposes of 
this argument that the equation of state is approximated by a 
7-law, i.e. P = (7— l)pe. The first and second assumptions 
completely determine the upstream flow as a function of ra- 
dius. Because the upstream flow is pressureless, we use the 
strong shock limit to determine the shock compression ratio, 
i.e. Pd/pu ~ (7+l)/(7 — 1)- Under these assumptions, the full 
Rankine-Hugoniot jump conditions reduce to a single expres- 
sion: 

Pd = p u J u (l-j^j (18) 

where u (d) denotes upstream (downstream) state variables. 
Since the term in parentheses is ~1, the shock conditions re- 
duce to an expression that demands a balance between the 
upstream ram pressure and the downstream thermal pressure. 
In essence, the primary independent shock condition is the 
momentum jump condition (eq.[To*ll. Form here on, we focus 
on the momentum jump condition with zero pressure on the 
upstream side, i.e. 

Pd+prt (19) 

Including turbulent pressure, the momentum jump condi- 
tion becomes 

Pel + PdV 2 d + PdRrr ~ PuV 2 u , (20) 

where the velocities, vj and v„, are background velocities. 
Hence, the new shock position is located where the post-shock 
thermal, ram, and turbulent ram pressures balance the pre- 
shock ram pressure. The addition of the turbulent ram pres- 
sure likely results in larger shock radii. 

Equation < [20b by itself is not enough to determine the shock 
position. One must also specify the pre-shock and post-shock 
profiles, and it is the intersection of these profiles that deter- 
mines the shock radius. The pre-shock ram pressure is given 
by free-fall assumptions, resulting in a fixed, relatively shal- 
low profile (e.g. p u v\ oc r" 5 ' 2 ). The post-shock region is in 
sonic contact and in rough hydrostatic equilibrium, so the 
postshock pressure depends upon physics (such as cooling) 
of the entire postshock region. Fortunately, though, the post- 
shock pressure profile can be expressed by a simple power- 
law (e.g. P ex r" (3_4) ), where the normalization depends upon 
the details of cooling, etc.. In section |H we fit power-laws to 
the pre-shock and post-shock profiles of 3D simulations and 
use eq. ( f2Qb to predict the average shock radius with and with- 
out turbulent ram pressure. We calculate the correct average 
shock radius only if we include the turbulent ram pressure. 
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Figure 2. Turbulent entropy luminosity , L s = 4nr 2 (F s ) vs. radius at 250 
ms after bounce for three driving neutrino luminosities (L„ = 2.1, 2.23, and 
2.3xl0 52 erg/s). In general, these profiles are consistent with a neutrino- 
driven convection hypothesis. L s is positive in the gain region where buoy- 
ancy actively drives convection, and it is negative where stabilizing entropy 
gradients cause buoyant deceleration. As is expected for neutrino-driven con- 
vection, the magnitude of the turbulent luminosity monotonically increases 
with the driving neutrino luminosity. 
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Figure 3. Reynolds stress as a function of radius and driving neutrino lu- 
minosity at 250 ms after bounce. All three diagonal components are shown: 
solid lines correspond to R rr , dashed lines correspond to Rgg, dot-dashed 
lines correspond to R^. On average, R rr ~ Rgg globally, and a 

Rgg locally. This equipartition in kinetic energy between the radial and tan- 
gential components is a commonly observed feature in buoyancy-driven con- 
vection, and is a consequence of buoyant driving in the radial direction, redis- 
tribution to the tangential components, and turbulent dissipation among all of 
the components. As is expected for neutrino-driven convection, R n - increases 
with neutrino luminosity. 

4. RESULTS 

In this section, we present many ways in which 3D CCSN 
simulations are consistent with the hypothesis that neutrino- 
driven convection dominates the multi-dimensional motions. 
Figure|2]shows the turbulent entropy luminosity, L s = 4irr 2 F s , 
versus radius for 3 different driving neutrino luminosities at 
250 ms after bounce. The gross features of these profiles are 
consistent with buoyancy-driven convection. In particular, L s 
is positive where convection is actively driven by heating in 
the gain region, and L s is negative where convective plumes 
overshoot into the cooling region and are decelerated by sta- 
bilizing entropy gradients. Furthermore, the magnitude of L s 
monotonically increases with the driving neutrino luminosity, 
an expected result for neutrino-driven convection. 

Similarly, neutrino-driven convection explains the 
Reynolds stresses. In Figure [3] we plot the radial (R rr , 
solid line) and tangential components (Rgg and R^, dashed 
and dot-dashed lines) of the Reynolds stress vs. radius at 
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Figure 4. G lobal buoyant driving, Wj,, vs. Global turbulent dissipation, E^. The fact that buoyant driving equals turbulent dissipation as predicted by turbu- 
lence theory ( Amett et al. 2009; Garaud et al. 2010; Murphy & Meakin 201 1 ) is a strong indicator that the neutrino-driven convection dominates the aspherical, 
nonlinear flow. 



250 ms after bounce. Two characteristics of the profiles are 
consistent with buoyancy-driven convection. First, like L s , 
the strength of the turbulent stresses (mostly R rr ) increases 
monotonically with neutrino luminosity. Secondly, the radial 
component of the turbulent stress is approximately equal to 
the combined tangential components; i.e. R rr ~ R^+Rgg, 
This result is consistently seen in other numerical ex- 
periments wh ere turbulence is unambiguously driven by 
buoyancy (see lArnett et al.l d2009l) , and r e ferences therein). In 
analytic derivation s (lArnett et all 12009b iGaraud et all 12010b 
iMurphy & Meakinl 1201 ll) . this approximate equipartition 
arises because buoyancy acts first on the radial compo- 
nent, and then the turbulence is dissipated after energy is 
redistributed among the three components. 

A prediction of the turbulent kinetic energy equation 
(eq. fTTT i is that, globally, buoyant driving is balanced by tur- 
bulent dissipation (eq.fl2l. Figure Invalidates this prediction. 
In this plot, we show global buoyant driving, Wt, vs. global 
turbulent dissipation, E^, for several driving neutrino lumi- 
nosities. In summary, Figure E] verifies that W b oc / pR% 2 dV, 
that the constant of proportionality is the largest scale as pre- 
dicted by Kolmogorov, and that global buoyant driving bal- 
ances global turbulent dissipation. 

In section [3] we argue that the convective luminosity, TqL s , 
and turbulent dissipation, £V, scale linearly with the driving 
neutrino-power, L v r. Figures|5]and|6]compare these analytic 
scalings, eqs. ( fOlfl4l i. (dotted lines) with the results from the 
3D simulations (diamonds). Indeed, both the turbulent lumi- 
nosity (Figure [3]l and turbulent dissipation (Figure |6) are lin- 
early proportional to the driving neutrino power. To calculate 
the constant of proportionality, we merely report the ratio of 



the luminosities at the last points. Coincidentally, the con- 
stants of proportionality, a and (3, approximately sum to one. 
From this coincidence, one might be tempted to conclude that 
ToL s + Ek ~ L v t. However, this supposition would contradict 
global balance of driving and dissipation which is a firm theo- 
retical prediction. Rather, we suspect that this interesting co- 
incidence is merely a consequence of the approximate deriva- 
tions leading to these scalings. Other than the fact that these 
constants are of order unity, we put no real emphasis on their 
values. The main result is that turbulence scales linearly with 
neutrino power, as is expected in buoyant convection. 

Finally, we find that turbulent ram pressure explains in 
part the expansion of the shock radius. As we argue in sec- 
tion [3] the position of the shock radius is essentially deter- 
mined where the pre-shock ram-pressure profile and the post- 
shock pressure profiles satisfy the momentum jump condition 
(eq. 1201 . In essence, the shock radius is established where 
the ram pressure of the upstream flow is balanced by the ther- 
mal, ram, and turbulent pressure of the downstream flow (see 
eqs. ( fl9ll20l i and the surrounding text in section |3). Figure 
[7] shows for one representative neutrino luminosity fits to the 
upstream ram pressure (dotted line) and downstream thermal, 
ram, and turbulent pressures (solid line), as a function of ra- 
dius. Using these pressure profiles, we calculate the shock to 
be located where the upstream and downstream fits cross. To 
identify where the shock would be without the turbulent pres- 
sure, we construct similar models for the downstream flow 
that exclude the turbulent pressure. For comparison we show 
the actual average shock radius ((R s )) from a 3D simulation 
(dot-dot-dot-dashed line). Without the turbulent ram pressure, 
the calculated shock radius is ~ 40 km short of the actual av- 
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Figure 5. Convective luminosity (TqL s ) vs. the driving neutrino power, L Ve r. The symbols show the maximum value of TqL s (restricted to the gain region) for 
three 3D simulations, all at 250 ms after bounce. For neutrino-driven convection, we analytically expect this convective luminosity to be linearly proportional 
to the driving neutrino power. For comparison, the dashed line shows thi s lin ear expectation. Our analytic calculation does not determine the constant of 
proportionality, so using the 3D simulations, we find that a ~ 0.7. See eq. 1 1 3 1 and the associated text for the derivation of the analytic scalings with neutrino 
luminosity. 
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Figure 6. Turbulent dissipation rate, E%, vs. the driving neutrino power, 
L„ f T. Similar to Figure[5] the diamonds show the 3D simulation results and 
the dash ed line shows the expected linear scaling with drivin g neu trino power, 
eq. U4t . More approximations are involved in deriving eq. ( I14K so the scal- 
ing in this plot is less robust than the scaling in Figure [3] Nonetheless, this 
plot shows that turbulent dissipation is roughly proportional to the driving 
neutrino power, as is expected for neutrino-driven convection. 

erage shock radius, but when we include the turbulent ram 
pressure, we accurately calculate the shock radius. 

Figure [8] shows the resulting shock locations as a function 
of driving neutrino luminosity. This plot shows the mod- 
eled shock radii, including turbulent ram pressure (solid line) 
and the modeled shock radii excluding turbulent ram pres- 
sure (dashed line). For comparison, we show the calcu- 
lated minimum (triangles), average (diamonds), and maxi- 
mum (squares) shock radii for 3D simulations all at 250 ms 



10" 



E 1° 2 



^ 10" 



10* 



10" 



(P+/3v 2 +pR rr , dilil 
P: 3D sim 




- R 



100 150 



200 250 300 
Radius [km] 



350 400 



Figure 7. Upstream (u) a nd do wnstream (d) momentum terms in the momen- 
tum shock condition, eq. )20t , and the resulting shock position. In essence, 
the average shock radius is established where the ram pressure of the up- 
stream flow (dotted line) balances the thermal, ram, and turbulent ram pres- 
sure of the downstream flow (solid line). The solid and dotted lines are fits to 
their respective flows. The dashed line shows the average pressure from the 
3D simulations. For comparison, we show the average shock radius ((S,)) 
of the 3D simulation (dot-dot-dot-dashed line). Neglecting the turbulent ram 
pressure, one would predict a smaller shock radius. The predicted and actual 
average shock radius is ~40 km larger than it otherwise would be without 
the turbulent ram pressure. In Figure [8] we use these curves to calculate 
the location of the average shock radius for several neutrino luminosities and 
show that it is a monotonically increasing function of L„ f and that turbulent 
pressure pushes the shock out to larger radii. 



after bounce. Including turbulent pressure in the post-shock 
profile predicts shock radii that agree with the measured av- 
erage shock radius. Excluding the turbulent pressure under- 
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predicts the average shock radius. On the other hand, exclud- 
ing the turbulent pressure gives shock radius predictions that 
are consistent with the minimum shock radii. This suggests 
that the minimum shock radii occur at places and times where 
the fluctuating turbulent motions are instantaneously negligi- 
ble. On average though, the turbulent motions are not negli- 
gible and influence the average shock radius. 

Showing that turbulent ram pressure causes expansion of 
the shock radius does not by itself prove that buoyant-driven 
convection is responsible for the expansion. Any instabil- 
ity that leads to turbulence would give a similar prediction. 
However, the dependency of the shock radii and the turbu- 
lent pressure on L Ur does strongly suggest the prominence of 
neutrino-driven convection. 

5. CONCLUSIONS AND DISCUSSION 

We have identified in section |4] five ways in which the as- 
pherical, nonlinear flow of 3D CCSN simulations is consis- 
tent with buoyancy-driven convection theory. For one, the 
turbulent luminosity is positive in the gain region where buoy- 
ancy actively drives convection, and the turbulent luminosity 
is negative in the stably-stratified region where buoyancy de- 
celerates convective plumes. Secondly, the radial component 
of the Reynolds stress is in rough equipartition with the tan- 
gential components, i.e. R n ~ Rgg+R^; this result is ob- 
served in other contexts of buoyancy-driven convection and 
is expected when convectio n is driven radially, but dissipated 
amon g all the components (|Arnett et al.1 12009b iGaraud et al.l 
12010b IMurphy & Meakinll201 ll) . Thirdly, we find that global 
turbulent dissipation is balanced by global buoyant driving. 
Fourthly, both the turbulent luminosity and turbulent dissipa- 
tion scale linearly with the driving neutrino power, a result we 
derive from buoyant convection theory. Finally, we calculate a 
consistent solution for the average shock radius only if we in- 
clude the turbulent ram pressure in the shock jump conditions. 
In essence, these results are consistent with the hypothesis that 
neutrino-driven convection is the dominant multi-dimensional 
instability in our 3D CCSN simulations. 

These results also suggest that there is no need to invoke any 
other multi-dimensional instabilities, in particular the SASI, 
to explain the aspherical, nonlinear motions. Of course, our 
results do not prove that the SASI is absent. Rather, they 
suggest that the nonlinear motions are merely consistent with 
buoyant convection and that if the SASI is present it mim- 
ics buoyant convection, or is subdominant. Numerous an- 
alytic and numerical studies have shown that if the condi- 
tions are right, a nonlinear SASI arises. However, these stud- 
ies were performed largely in the absence of neutrino heat- 
ing. A few did include neutrinos, but as far as we can tell, 
neutrino-driven co nvection eventually dominates the flow of 
these simulations dFoglizzo et alj 12006b IScheck et all 12008b 
iFernandez & Thompsonll2009h . It appears that in attempting 
to isolate the SASI mechanism, researchers were suppress- 
ing the dominant nonlinear instability in core-collapse simu- 
lations. 

Minimally, there are two hypotheses for the dominant 
multi-dimensional instability: the SASI and buoyancy-driven 
convection. Neither hypothesis can be proven; they can only 
be shown to be consistent or disproven. Because a self- 
cons istent theory exists for n onlinear, buoyant convection 
(e.g. IMurphy & Meakinll201 lb . we can show that buoyancy- 
driven convection is consistent with post-shock turbulence. 
However, there is as yet no self-consistent theory for a non- 
linear SASI mechanism, so we can not test the SASI at this 



time. In fact, it is entirely possible that some aspects of 
nonlinear convection theory would be mirrored in a nonlin- 
ear SASI theory. For example, in both cases, turbulence is 
most likely dissipated in accord with Kolmogorov's hypoth- 
esis. However, we have shown that turbulent dissipation is 
balanced by buoyant driving, which seems an unlikely pre- 
diction of a nonlinear SASI theory. In any case, it is clear 
that a nonlinear theory for the SASI must be developed be- 
fore we can definitively claim that the SASI is subdominant. 
Furthermore, even though our approximations are designed to 
closely mimic more self-consistent simulations, a robust con- 
clusion on the importance of convection must wait for full 3D 
neutrino-transport hydrodynamic simulations. 

In the meantime, we have shown that a theory for nonlin- 
ear, buoyant convection is consistent with the post-shock tur- 
bulence of 3D simulations that include neutrino heating. Not 
surprisingly, we find that neutrino-driven convection accom- 
panies neutrino-driven explosions. 
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